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O ! Abstract 

^ : 

^ ■ We compute the otluf perturbative QCD contribution to semileptonic B decay, in- 

^ \ eluding the finite mass of the charm quark. This result provides an estimate of the size of 

. the two-loop correction, which is found to be about 50% of the one-loop correction. We 
OO ■ 

' use these results to set the scale for the one-loop correction using the scheme of Brod- 

o : 

^ . sky, Lepage and Mackenzie and find a BLM scale of /Ublm = 0.13 m;,, when the inclusive 

Q\ • semileptonic rate is expressed in terms of the h and c quark pole masses and the MS strong 

\ coupling. The two loop correction lies roughly midway between that obtained at rUc = 
P-i! 

p . and that obtained in the Shifman-Voloshin limit m^, rric ^ — rric ^ Aqcd while the 

(D : 
^ ■ 

•1—1 , 
X- 

^3: 



corresponding BLM scale is somewhat closer to that obtained in the former case. 



UTPT 94-27 
CMU-HEP 94-32 
CALT-68-1955 

DOE-ER/40682-86 October 1994 



1 



1. Introduction 



Inclusive semileptonic B decays provide a method for determining the magnitude of 
the element of the Cabibbo-Kobayashi-Maskawa matrix Vcb- In the limit where the b 
quark mass is much larger than the QCD scale the B meson decay rate is equal to the b 
quark decay rate [0. Corrections to this first arise at order [h-Qco/'n^bY these nonper- 
turbative corrections may be written in terms of the matrix elements {B\b{iD)'^b\B) 
and {B\biga ^jjG^^b\B) . The order as corrections to the decay rate, including the complete 
c quark mass dependence, have been calculated in |^. 

No complete calculation of the two-loop perturbative corrections to inclusive heavy 
quark decay has been performed. However, the 0{a1nf) correction to the inclusive semilep- 
tonic decay rate to a massless quark has been recently calculated 0. This correction is 
of interest in its own right for two reasons. First, if one adopts the viewpoint of Brodsky, 
Lepage and Mackenzie [[7|, one can use this result to set the appropriate scale for ctg in 
the leading term. For b Xu semileptonic decay, this leads to the surprisingly low scale 
A*BLM = 0.07 mfo 0. Second, even if one does not adopt this viewpoint, we argue that 
this calculation provides a good estimate of the size of the two-loop calculation. This is 
not because n/ is large, but rather because the QCD beta function /3 = 11 — 27i//3is 
large, so the vacuum polarization graphs which contribute to this term are be expected to 
dominate the two-loop result. Empirically, this is certainly true for R{e~^e~ hadrons). 



r(r — > z/^- + hadrons) and the two-loop relation between the pole mass and the MS mass 
of a heavy quark: 



R{e+e~ hadrons) = 3 ( ^ Q,^ 



l + + (0.17/3 + 0.08) 

TV \ TV / 



r(r Uj- + hadrons) as{mr) i \ ( «s(w. 



_ 1 + ^^ + (0.57/3 + 0.08)^^ +... 
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where ols is the MS strong coupling. In each of these cases, the 0(ag/3) term provides an 
excellent approximation to the full two-loop correction. Using the term proportional to /3 
as an estimate of the size of the two-loop QCD corrections for b —>■ Xu semileptonic decay 



and taking n/ = 3^ gives 



,2 '^F"''b 

' 1927r3 



1- ( ^i^j [2.41] -f^i^]' [28.7] + ... 



(1.2) 



where we have used as(mb) ~ 0.2. The two-loop term is clearly significant. 

For charmed hadrons the situation is even worse. Using as(mc) = 0.29, we find 

r(c XaeiT^) = \V,d\^-£-f[l - 0.22 - 0.25 + ...] (1.3) 

and it does not appear that the perturbation series is controlled at all. This is reflected in 
the small size of the corresponding BLM scale /Ublm ~ 100 MeV. 

The results of |j|] raise questions about the applicability of perturbative QCD to 
describe the semileptonic decay of charmed hadrons, and also suggest that two-loop cor- 
rections are important for a reliable extraction of Vbc from inclusive b —>■ Xc decays. In this 
paper we extend these results to include the effects of the finite c quark mass in 6 — > Xc 
decays. We find that, for mc/rrib = 0.3, the BLM scale is raised to /Ublm = 0.13mb, and 
the corresponding perturbation series is 



n2 5 

T{b^XceV;) = \V,c\^-I-±[0.52] 



1 - ( ^i^) [1.67]- (^i^)' [15.1] + ... 



(1.4) 

The finite charm quark mass thus reduces the size of the higher-order corrections, although 
they are still significant. 



In Eqs. (p..2| ) - (|L4D , rrib and rric are the pole masses of the b and c quark. If MS masses 
are used in ( |1.2| ) and ( prB| ) the convergence is slightly improved [^. 



^ We have not taken into account the c quark mass for virtual c quarks. The true result will be 
somewhere between n/ = 3 and n/ = 4, and we have chosen to take the number of light flavours 
to be three for both c and b decays. 
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2. Calculation 



The semileptonic decay rate of a 6-quark is given by 



n / 3 \ 4 \mb 



Tib ^ X^eiT-,) = \Vbc\' f( — 



n J \ \mb / \mb 

' ' (2. 

where /9 = 11 — |?t./, is the MS' strong coupling, f{x) = (1— x^)(l— 8x^+x^)— 24x'* log(x) 
and mb,c are the heavy quark pole masses. The function di{x) has been computed in f^, 
and 5i(0.3) = —1.11. In this work we compute the function xpi^)- 

The second order corrections to the semileptonic decay rate have been shown to 
be calculable from the first order corrections when a finite gluon mass is included. If the 
first order correction with a finite gluon mass rUg is denoted by T^^\mg) then the second 
order correction proportional to the perturbative QCD beta function is 

where a^^^ is the strong coupling constant evaluated in the "V-scheme" of [0 and is related 
to the MS coupling constant ag by 

a(/Hp)=as{p) + m^f3 + .... (2.3) 

Both virtual and bremstrahlung graphs contribute to the first order correction. The 
virtual graphs were computed analytically, with the final integrations of the lepton q"^ 
and mgiuon performed numerically, while the phase space integrals for the bremstrahlung 
graphs were performed numerically. 

The leading two-loop correction xpij^c/'n^b) is plotted in fig. According to the BLM 
scheme [0, this term is absorbed into the first order correction by redefining the scale at 
which the first order result is evaluated to n = /Ublm, where 



AtBLM = mb exp 



7r2 - ^ + Si{mc/mb) 



(2.4) 



The BLM scale is plotted as a function of rric/mb in fig. |^. At = we reproduce the 
result //blm(O) = 0.07 m^, while for rric/mb = 0.3 we find /iBLM = 0.13 mf,. This is 



somewhat higher than the scale found for the case of massless final state quarks but is still 
considerably smaller than m^. 

A non-trivial check on our numerical results is obtained in the Shifman-Voloshin limit 
mh, rric ^ mi, — rric ^ Aqcd 0- In this limit the final hadronic states are dominated by 
the D and D* mesons (which are degenerate to leading order in l/rric)- They are produced 
almost at rest, and the rates for these exclusive processes can be calculated in the heavy 
quark effective theory with no free parameters. This gives 



Gpjmb - rric 
607r3 



Owl' + 3|r7A|'] 



nib 



(2.5) 



where the coefficients rjA and rjv arise in matching from the QCD currents to the heavy 
quark currents 

cYh^TlvKl^hb + ... 

(2.6) 

C7^75^ VAhcl^l->hb + ... , 
and in this limit the heavy c and h quarks have the same 4-velocity. The matching coeffi- 
cients have the perturbative expansion 



r]v 



1 + (f){mc/mb) + 

6 TV 



OLs{mb) 



TV 



2 r 



1 



[mc/mb)P + 



+ 



1 + 0(mc/mb) - 2 + 

6 n 



as{mb) 



IT 



-^(j){mjmb) - ^ ) /3 + ... 



+ ... . 

(2.7) 

where the leading term was calculated in Refs. and [jl^ and the 0{a'^P) terms were 
calculated in Ref. []TT[|. The function (f){z) is defined by 

1 + z 



1 - z 



log z — 6 



21 



and so ^(1) = 0. Furthermore, at rric — "^6 the vector current is not renormalized, and 
Eq. ( |2.7D simplies to 

= 1 



1 



2 as{mb) 7 f as{mb) 



n 



36 



TV 



+ ...] + 



(2.9) 



which gives 



TiB-.D^D*el7c)=^'^^'^'-'^^^' 



+ 



157r3 

{rrib - rric)^ 
rrib 
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The one loop term is in agreement with the results of Jezabek and Kuhn 0, while the 
two-loop term is in agreement with the limiting value of our inclusive calculation. This 
gives a BLM scale in the SV limit of hblm = 0.56 mb. 



3. Discussion and Conclusions 

The BLM scale is one measure of the size of higher order corrections. However, as we 
argued in the introduction, we can also use our results as an estimate of the size of the 
two-loop corrections with /x = nib- For mc/rrih = 0.3, the inclusive rate at leading order in 

l/nib is 



r<2 5 

r(6^X,ez7^) = |y,,|2^g|[0.52] 



as{mb) \ „_i / OLs{mb 



1 - ( ^^^^ 1 [1.67] - ( ^^^^ ) [15.1] + ... 



2 



r<2 5 

= lHcp-rT^[0.52] [1 - 0.11 - 0.06 + ...] 

(3.1) 

where the factor of 0.52 arises from the lowest order phase space factor /(0.3) and again we 
have taken ^^(mb) = 0.20. The one and two loop corrections relative to the leading term 
are plotted as functions of rric/mb in fig. |^. We can compare this with our previous results 
for b —>■ Xu semileptonic decay, Eq. (O). The finite charm quark mass somewhat reduces 
the size of the two-loop term, but it is still significant, about 50% of the 0{as) correction. 
The two loop correction lies roughly midway between that obtained at rric = and that 
obtained in the Shifman-Voloshin limit while the corresponding BLM scale is somewhat 
closer to that obtained in the former case. 

We have expressed the b Xc decay rate in terms of the b and c quark pole masses. 
This is convenient since the difference in these masses is determined by the meson masses 
(up to corrections of 0{l/mc,b))'- 

nib - rric = ruB - ruD + O {l/'mb,c) (3.2) 

and the function f{mc/mb) is less sensitive to the sum of the quark masses than the 
difference. 

For charm decays, including the finite mass of the s quark has a negligible effect on 
the magnitude of the two-loop corrections. Because the s quark is so light its pole mass 
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is not a useful quantity. However, for illustrative purposes we have plotted the size of the 
one and two loop corrections relative to the leading term as functions of mg/mc in fig. ^ 
using as(mc) = 0.29. Note that for small nis, the two loop correction is larger than the 
one loop correction. For ms/rric = 0.13 the inclusive semileptonic decay rate is 



r(c — > XgeUe) 



^^'^ 1927r3 

^ r.! 



[0.87] 



1 - 



^) [2.08] - f ^ 



[22.7] + 



(3.3) 



[0.87] [1-0.19-0.19 + 



1927r3 

The perturbation series is still clearly uncontrolled. 
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Figure Captions 

Fig. 1. The function X/si^c/^b) from Eq. plotted vs. rric/mi,. 

Fig. 2. Tlie BLM scale as a function of rric/mb- 

Fig. 3. The 0{as) (dashed line) and 0{a'^(3) (solid line) corrections to the expansion of 
r(6 XcCVe) relative to the leading term, plotted as functions of mc/rrih- We 
are using aairrib) = 0.20 and n/ = 3. 

Fig. 4. The 0{as) (dashed line) and 0{a'^(3) (solid line) corrections to the expansion of 
r(c — s> XgeVe) relative to the leading term, plotted as functions of ms/rric. We 
are using asirric) = 0.29 and n/ = 3. 
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